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Abstract. A quasiparticle description of pure glue QCD thermodynamics at T < Tc is proposed and 
compared to recent lattice data. Given that a gas of glueballs with constant mass cannot quantitatively 
reproduce the early stages of the deconfinement phase transition, the problem is to identify a relevant 
mechanism leading to the observed sudden increase of the pressure, trace anomaly, etc. near the critical 
temperature. It is shown that the strong decrease of the gluon condensate at T < Tc combined with the 
increasing thermal width of the lightest glueballs might be the trigger of the phase transition. 
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1 Introduction 

The phenomenology related to the deconfinement phase 
transition from hadronic matter to quark-gluon plasma 
(QGP) at high enough temperatures or densities is a very 
active field of research since the beginnings of QCD, with 
pioneering works such as Refs. [BIl]. Much effort is de- 
voted to understand the QGP on the experimental side: 
The QCD matter is or will be studied in heavy-ion colli- 
sions at RHIC, SPS, FAIR, and the LHC (see Refs. HE] 
for more information). On the theoretical side, the QGP is 
also a challenging topic, which has deserved lots of stud- 
ies within several frameworks. In particular, phenomeno- 
logical quasiparticle models and more fundamental ap- 
proaches like perturbative calculations or AdS / CFT dual- 
ity have proved to be successful both in reproducing lattice 
data and in providing a reliable description of experimen- 
tally observed results at RHIC, see e.g. Refs. [5l|6l[7] for 
some recent results and Refs. [S] for some reviews on the 
topic. 

Lattice QCD has a particular status since it is, in prin- 
ciple, the most powerful technique to extract nonpertur- 
bative informations from QCD. When experimental data 
are lacking, lattice data are then often used to fit other 
model's parameters. A crucial observable in QGP physics 
that can be computed in lattice QCD is the QCD equa- 
tion of state (EoS). In particular, the EoS of pure glue 
SU(3) QCD has been computed on the lattice more than 
a decade ago [9j, while more detailed data have recently 
been obtained in Ref. [7] with gauge groups ranging from 
SU(3) to SU(8). Many other calculations of the QCD EoS 
have also been performed including quarks flavors, at zero 
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chemical potential or not. The interested reader will find 
more references in e.g. the review (10) . 

Among the various existing phenomenological approaches, 
quasiparticle models rely on the assumption that the QGP 
can be seen as a gas of deconfined quarks and gluons. Since 
pioneering works , they have been shown to be able to 
reproduce the various EoS computed in lattice QCD at 
T > Tc. It is thus tempting to apply such models below 
Tc- In that case however, the QCD matter should rather 
be modeled by a hadron gas since QCD it is in its confin- 
ing phase. We focus here on the pure glue EoS below Tc, 
which already captures the essential physical features of 
the full QGP without involving extra technical difficulties 
due to the presence of quarks. 



2 Glueball gas model 

The pure glue hadronic matter is expected to consist in 
glueballs and, since the scattering amplitudes between 
glueballs scales in l/iV^ (instead of for mesons), a 

noninteracting Bose gas may be assumed in a first ap- 
proach. Following Bose-Einstein statistics, only the lowest- 
lying glueballs should bring a significant contribution to 
the EoS since the one of a glueball of mass nig and spin Jg 
is roughly proportional to {2Jg -\- l)e~™s/"^. However, by 
using typical values for the low-lying glueball masses |12| 
in a hadron gas model, one fails to reproduce the strong 
increase of the thermodynamical variables 
ready pointed out in Refs. P^fTl] . To our knowledge, the 
only proposal leading to a model in agreement with the 
lattice EoS is the one of Ref. [13] , where a high- lying glue- 
ball spectrum of Hagedorn-type is assumed. The expo- 
nential growth of the number of states with mass mg then 
compensates the statistical suppression and brings a large 
enough contribution near Tc. It is also worth mentioning 
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the recent Ref. [TB], where it is shown that a Hagedorn- 
type model of the QGP is able to reproduce the trace 
anomaly computed in lattice QCD while ensuring a low 
viscosity/entropy ratio. 

The Hagedorn spectrum relies on a string-theoretical 
picture of the hadron spectrum. In particular, following 
the ideas of Ref. [13], one has to assume that glueballs 
are excitations of a closed string, which is not so obvious, 
partly in view of the many models reproducing the lattice 
spectrum at zero temperature by assuming totally differ- 
ent frameworks [13j. That is why it might be of interest 
to see whether an alternative way of understanding the 
early stages of the pure glue phase transition (T < T^) 
can be found or not. This is the purpose of the present 
work, where it is chosen not to assume a Hagedorn glue- 
ball spectrum. Therefore, only the lightest glueballs, that 
all models find to be the scalar and tensor ones in pure 
glue QCD [IIT, are assumed to bring a significant contri- 
bution to the EoS. 

As it will be shown, two ingredients are needed. First, 
the nontrivial contribution of the gluon condensate to the 
trace anomaly, especially its brutal reduction near the 
phase transition. Remark that the strong link between the 
gluon condensate and the QCD phase transition have been 
previously outlined for example in Refs. [T7]. Second, a 
significant reduction of the glueball masses near Tc. Since 
that point is the most counterintuitive one at first sight, 
let us analyze it first. 
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Fig. 1. Scalar and tensor glueball masses computed through 
a pole- mass fit in Ref. [19] (symbols), compared to a fit of the 
form (dJ with (in units of Tc) To++ = 0.596, m°++ = 5.547, 
bo++ = 4.230, T^++ = 0.755, ml++ = 8.113, and 63++ = 7.152. 
A standard value Tc = 265 MeV is used. 

are well described below Tc by the form ([Ta|) with 

TOg(r) = m°, 

r<,(r) = o T<Tg 

= bgiT~ Tg) Tg<T< Tc . (lb) 



3 Glueball mass reduction 

The behavior of glueballs at finite temperature has been 
first investigated on the lattice in Refs. [T51IT5] . The basic 
conclusion of those works is twofold, following the proce- 
dure used to analyze the results. On one hand, if the tem- 
poral glueball correlator is fitted assuming glueball states 
with zero width, then the 0++ and 2++ pole masses are 
found to significantly decrease for T < Tc- On the other 
hand, a Breit-Wigner fit can be applied to that corre- 
lator; then the glueball masses are found to be roughly 
constant from T = to Tc, with a thermal width which is 
nearly zero at low temperature but then increases rather 
linearly near Tc- That conclusion is quite intuitive and 
might be a general feature of glueballs near Tc- The glue- 
balls should progressively be "dissolved" in the medium 
when approaching the deconfinement temperature, and 
their width should increase accordingly. Notice that such a 
constant mass is expected from calculations relying on ef- 
fective low-energy Lagrangians |20] . It is moreover shown 
in Ref. 19 that the pole mass, nig{T), and the Breit- 
Wigner mass, mg{T), and thermal width, rg{T), are linked 
as follows: 

TUg (T) « mg (T) -2T+ ^^T^-TgiTy. (la) 
As shown in Fig.[Tl the pole masses computed in Ref. [TO] 



Such an ansatz qualitatively encodes the results of the 
Breit-Wigner fit performed in Ref. [TH]. Remark that the 
thermal broadening of the glueballs generates a pole-mass 
reduction; interestingly such a glueball mass reduction 
near the phase transition has also been predicted by com- 
putations within the dual Ginzburg-Landau theory |21) . 
We point out that the value Tc = 265 MeV will be as- 
sumed in the rest of this paper; it lies in the typical range 
260 - 280 MeV found in pure glue QCD. 



4 Glueball gas pressure 

As said in Sec.j^l we assume the interactions between glue- 
balls to be negligible in a first approximation. The pres- 
sure of an noninteracting gas of glueballs with mass rUg 
and spin Jg simply corresponds to that of an ideal Bose- 
Einstein gas [35] 

p^^Jllll^ dke\n{l-e-V^^'^). (2) 

Such a framework is the most straightforward way to 
deal with a boson gas. However, Eq. ^ shows that only 
the glueball masses appear in the computation of the EoS, 
and we have seen in the previous section that the in- 
creasing thermal width should also be taken into account. 
Rigorously, one would thus need a statistical formalism 
not only generalizing the Bose-Einstein distribution to the 
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case of unstable particles, as done first in [21], but also ex- 
plicitly including the interactions between glueballs, since 
an increasing width should enhance those interactions. 
Such a formalism is out of the scope of the present paper, 
which should be seen as a first step towards an alternative 
way of understanding the pure glue EoS below Tc- That 
is why two assumptions will be made in the following. 
First, we still choose to neglect the interactions between 
glueballs. As discussed in the introduction, large- iVc argu- 
ments might be relevant to justify that point. Second, we 
consider that the nonzero-width effects can be absorbed 
into a redefinition of the glueball masses. Indeed, we have 
shown in the previous section that a glueball with a con- 
stant mass and an increasing thermal decay width can 
be described effectively as a glueball with a zero width 
and a decreasing pole mass. That is why the pole mass 
rag{T) will be used hereafter. Notice that the second ap- 
proximation is expected to be less and less accurate when 
approaching Tc from below, since the ratio rg(T)/mg{T) 
is an increasing function of T. From [TH] one can compute 
that ro++{Tc)/mo++ « 0.34 and r2++{Tc)/m2++ « 0.25, 
meaning that the glueball decay widths become significant 
with respect to their masses at the phase transition. Nev- 
ertheless one can reasonably think that our approach will 
be able to predict the global trend of the EoS near Tc. 

Starting from Eqs. ^ in which — mg(T) given 
by (|lap . the pressure of the QCD matter below Tc should 
then be given by p = J2gPg^ the sum running over all 
the glueball states. But, since a Hagedorn spectrum is 
not assumed, only the lowest-lying glueballs will signif- 
icantly contribute because the statistical suppression in 
{2Jg + l)e~'"s/^ is not balanced by the exponentially ris- 
ing number of states with respect to nig. Moreover, the 
higher mass glueballs would contribute if rUg/T was of or- 
der unity, but rUg/Tc » 1 even for the 0++ glueball. We 
thus take 

p«Po+++P2++- (3) 

Results obtained from this last equation are shown in 
Fig. [5] and compared to the lattice data of Ref. (Tj . The 
conclusions are the following. First, using a glueball gas 
with constant masses nig ~ ni^ fails to reproduce the ob- 
served increase of pressure near Tc, even by taking lighter 
values than the fitted ones, given in Fig. [T] Second, using 
the temperature-dependent masses, nig = nig{T), fitted 
from Ref. [T9l , greatly improves the agreement with lattice 
QCD and is a first argument in favor of the scenario pro- 
posed here. Remark that p{Tc) is underestimated; we will 
show in the next section that the contribution coming from 
the gluon condensate is able to cure this problem. Another 
possibility could be that the higher-lying glueballs bring a 
significant contribution to the pressure. Provided that the 
mass reduction mechanism does not cause those glueballs 
to become lighter than the tensor one, it can be computed 
that the contribution shifts the pressure of less than 
7% (the ^ state is the closest in mass of the 2++), while 
the heavier states contribution is even smaller. Thus only 
the scalar and tensor glueballs may be considered in a first 
approximation. It is worth mentioning the recent work [2 3) 
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Fig. 2. Pressure computed in pure glue SU(3) lattice QCD 
and normalized to the Stefan-Boltzmann factor k,sb ~ 45/87r^, 
taken from Ref. J (full circles) . The lattice data are compared 
to Eq. (131 using either nig = mPg (dashed line), nig = mg{T) 
given by Eq. ([1} (dotted line) , and the full model in which the 
gluon condensate contribution ((8)1 is added (solid line). 



suggesting that, indeed, the + glueball remains heavier 
than the tensor one. 



5 Gluon condensate contribution 

The next step is now the computation of the trace anomaly, 
a priori straightforwardly defined from the pressure Q by 

^ = T'9t{^). (4) 

A look at Fig. [3] clearly shows that our model with mg{T), 
although satisfactorily reproducing the lattice pressure, 
severely underestimates the trace anomaly near the phase 
transition. This seemingly paradoxical situation can be 
clarified as follows: When speaking of the trace anomaly, 
the nontrivial role of the gluon condensate has to be taken 
into account. That argument has already been proposed 
when studying the QCD matter at T > 1.2 Tc [3S]. Let us 
now apply it below Tc. 

It is known that the gluon condensate at temperature 
T, defined by (G^)^ = - (f G;i^G(^''(T)^, contributes to 

the QCD trace anomaly as = (G^) - (G^)^ 
Thus the total trace anomaly. A, should rather be 

A = A + Ag2. (5) 

The gluon condensate can moreover be written as the 
sum of a magnetic and an electric part, i.e. (CP^^ — 
(Gc)^-I-(G^)^ in euclidean space. It appears from lattice 
QCD simulations that (G^)^ « (G^)^ on one hand, and 
that {Gl)^ is such {Gl)^ « {G^)^ /2 but then falls very 
quickly near Tc to reach a zero value just after the phase 
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Fig. 3. Same as Fig. [5] for the trace anomaly. The dashed line 
is the glueball contribution (|4]), while the solid line comes from 
Eq. (O, where Aqi = Aq2 following Eq. Q. 



transition [22] ■ Consequently, one expects [25] 



Ag2 = 



[l-Ce(T)], where Ce(T) = 



(6) 



can be known from lattice computations [2 7) . 

Since Ce{T) is only known at a few temperatures from 
the lattice (see Fig. S]), it would not be convenient to use 
the available data as an input in our model. Rather, it is 
better to compute first the values of Aq2 that fit the trace 
anomaly obtained in Ref. [7| and then to check whether 
they agree or not with the results of Ref. ^7} . The values 
needed to fit the trace anomaly computed in Ref. [7] is 
very accurately fitted by the form 



1 - 



T 



1 



fl + e(T/Te-/3)//4^ 



T > 0.9 Tc, 
r<0.9rc, (7a) 



with 



/i = 6.171, /2 - 0.637, h = 1-013, U - 0.015. (7b) 

As shown in Fig. [31 the trace anomaly ([5]) computed 
with Ag2 = Aq2 fits the lattice data very well - it is log- 
ical since it has been designed for. The key observation is 
now that, in Fig. (H the fitted term © is in good agree- 
ment with that obtained by putting the values of Ce(T) 
available from the lattice study [37] in the theoretical es- 
timate ©. A typical value (G^)^ =0.030 GeV* has been 
used as a mean value of various results obtained so far in 
the literature [25] . Such an agreement between the needed 
value of the gluon condensate and the one theoretically ex- 
pected is a relevant check of the mechanism presented here 
describing the phase transition. 
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Fig. 4. Trace anomaly contribution Aq2 needed to fit 
the lattice data presented in Fig. [3] (solid line). The curve 
agrees with formula ^ in which c^iT) has been taken 
from the lattice study [57], together with the typical value 
(G^)^ =0.030 GeV* =6.08 T," (full circles). 



The coherence of our scenario requires the gluon con- 
densate contribution to the pressure to be computed from 
the thermodynamical relation Q and then to be added 
to the glueballs pressure p. One has 



Pg2 = T 



^ Ag^{x) 



dx. 



'0 ^ 

and the total pressure finally reads 



P = Po++ +P2++ +Pg^- 



(8) 



(9) 



It is readily observed in Fig. [2] that, near Tc, the total 
pressure (O reaches an excellent agreement with the lat- 
tice data. In particular, the pressure in T = Tc is no longer 
underestimated. 

Since in Ref. [7 , to which our model is compared, 
the energy density as well as the entropy density are ob- 
tained as linear combinations of the pressure and trace 
anomaly following standard thermodynamics, it is enough 
for our purpose to have considered p and A. We neverthe- 
less mention for completeness that the entropy density, 
s — OtP, has been independently computed on the lat- 
tice in Ref. [M] between 0.7 Tc and T^. The present model 
actually overestimates the dimensionless ratio s/T'^ ob- 
tained in this last work, especially near T^- We hope that 
more lattice data concerning the QCD EoS below Tc will 
be available in the future in order to refine our model and 
perform a more quantitative comparison with the different 
lattice results. 



6 Large Nc limit 

Some remarks can be done about the large- A'^c limit of our 
model, that has not been considered up to now. By defini- 
tion, KsB oa 1/{N^ — 1). Above Tc, the dominant contribu- 
tion to the equation of state should come from deconfined 
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Fig. 5. Total pressure @ normalized to N^T'^ computed at 
Nc — 3 (solid line) and at large Nc (dashed line). 

gluons, whose {N^ — 1) color degrees of freedom cancel the 
KSB factor, leading to globally constant thcrmodynamical 
properties with respect to Nc as observed in Ref. [7 . Be- 
low Tc however, all glueballs are in a color singlet and their 
mass scales as = 0{1). The glueballs contribution to 
the EoS normalized according to ksb should thus tend to 
zero at large Nc- But, since the gluon condensate scales 
as (Nc — 1), its thermodynamical contribution remains 
constant at large Nc. The trace anomaly and pressure at 
large Nc are then expected to read A w Aq2 and p w pq2 
respectively in our model. 

Following our scenario, the gluon condensate term is 
dominant in the trace anomaly but not in the pressure. 
Thus, we predict that the pressure below Tc will be strongly 
suppressed at large Nc as illustrated in Fig. [51 while the 
trace anomaly will decrease only slightly. The results of 
Ref. [2 suggest that the trace anomaly could be more 
and more peaked near Tc at large Nc, leading to a phase 
transition which is more of first-order type, in qualitative 
agreement with what is expected from the present ap- 
proach. Remark however that recently, a strong expansion 
of Yang-Mills theory at large Nc has been proposed [2£| , 
in which the pressure is found to scale as N^, with no 
Nc term as in our model. It can be hoped that future 
studies of glueball and glueball condensate properties at 
finite-temperature and large Nc will allow a more accurate 
validation of the ideas developed in this work. 



7 Conclusion 

A new way of understanding the pure glue QCD equa- 
tion of state below Tc has been proposed. The basic idea 
is that the pure glue hadronic matter consists in both a 
glueball gas and a gluon condensate part. It can be intu- 
itively expected that the thermal width of glueballs tends 
to increase when approaching the deconfinement phase 
transition. Such a behavior is effectively taken into ac- 
count in our formalism through a decrease of the glueball 
masses near Tc- Moreover, that effect has to be combined 



with the vanishing of the gluon condensate at the critical 
temperature. 

In order to illustrate qualitatively the proposed sce- 
nario, lattice data have been used as numerical inputs in 
the model: The glueball masses and the gluon condensate 
values have been taken from Refs. [15] and [27] respec- 
tively. Those data, when incorporated into computations, 
lead to an equation of state in good agreement with that 
computed in Ref. [7] , thus providing an a posteriori coher- 
ent interpretation of various independent existing results 
in finite-temperature QCD. 

Remark that, although the agreement reached by our 
model and the lattice data is excellent, the present ap- 
proach aims at being a first step towards a more com- 
plete model. Its weakness is mostly the way of including 
the glueball decay width, i.e. by redefining the glueball 
masses without modifying the form of the Bose-Einstein 
distribution, which can be questionable very close to Tc 
where the decay widths become significant with respect 
to the glueball masses. The elaboration of a more realistic 
approach is left for future works. 

Finally, at large Nc, we predict a strong reduction and 
a weaker decrease of the pressure and trace anomaly be- 
low Tc respectively. The present discussion suggests that 
new high-precision computations of glueballs and gluon 
condensate properties as well as of gluon plasma thermo- 
dynamics below Tc might be worth in order to check the 
general scenario developed here and the suggested large- 
Nc behavior in particular. 
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